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ABSTRACT 

We report the first detection of an VIII Lya absorption line associated 
with an overdense region in the intergalactic medium (IGM) along the sight- 
line towards PKS 2155-304 with the Chandra Low Energy Transmission Grat- 
ing Spectrometer (LETGS). The absorption line is detected at 4.5a level with 
cz « 16,600 km s^ 1 . At the same velocity Shull et al. (1998) detected a small 
group of spiral galaxies (with an overdensity of 5 ga i ~ 100) and low metallicity 
H I Lya clouds. We constrain the intragroup gas that gives rise to the O VIII 
Lya line to a baryon density in the range 1.0 x 10~ 5 < < 7.5 x 10~ 5 cm -3 
(50 < 5b < 350) and a temperature of 4 — 5 x 10 6 K, assuming 0.1 solar abun- 
dance. These estimates are in accordance with those of the warm/hot inter- 
galactic medium (WHIM) that are predicted from hydrodynamic simulations. 
Extrapolating from this single detection implies a large fraction of the "missing 
baryons" (~ 10%, or ~ 30-40% of the WHIM) are probed by the O VIII absorber. 

Subject headings: intergalactic medium — quasars: absorption lines — large- 
scale structure of universe — BL Lacertae objects: individual (PKS 2155-304) 

1. Introduction 

The cosmic baryon budget at low and high redshift indicates that a large fraction of 
baryons in the local universe have so far escaped detection (e.g., Fukugita, Hogan, & Peebles 
1998). While there is clear evidence that a significant fraction of these "missing baryons" 
(between 20-40% of total baryons) lie in photoionized, low-redshift Lya clouds (Penton, 
Shull, & Stocke 2000), the remainder could be located in intergalactic space with temper- 
atures of 10 5 — 10 7 K (warm-hot intergalactic medium, or WHIM). Resonant absorption 
from highly-ionized ions located in the WHIM gas has been predicted based on both ana- 
lytic studies of structure formation and evolution (Shapiro & Bahcall 1981; Aldcroft, Elvis, 
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McDowell, & Fiore 1994; Perna & Loeb 1998; Fang & Canizares 2000) and cosmic hydrody- 
namic simulations (Hellsten, Gnedin, & Miralda-Escude 1998; Cen & Ostriker 1999a; Dave 
et al. 2001; Fang, Bryan & Canizares 2002). Recent discovery of O VI absorption lines by 
the Hubble Space Telescope (HST) and the Far Ultraviolet Spectroscopic Explorer (FUSE) 
(Tripp & Savage 2000; Tripp, Savage & Jenkins 2000; Tripp et al. 2001) indicates that there 
may be a significant reservoir of baryons in O VI absorbers. While Li-like O VI probes 
about ~ 30 - 40% of the WHIM gas (Cen et al. 2001; Fang & Bryan 2001), the remaining 
~ 60 — 70% is hotter and can only be probed by ions with higher ionization potentials, such 
as H- and He-like Oxygen. 

The WHIM gas, predicted by hydrodynamic simulations, typically has an overdensity 
of 5b = 5 — 200 and is distributed in small groups of galaxies or in large scale filamentary 
structures that connect virialized halos. Given the expected physical conditions, current 
instruments such as Chandra and XMM should be capable of detecting resonance features 
from the WHIM gas (Fang, Bryan & Canizares 2002). A first attempt by Fang, Marshall, 
Bryan, & Canizares (2001) with the Chandra High Energy Transmission Grating Spectrome- 
ter (HETGS) yielded an upper limit of N(0 VIII) < 10 17 cm -2 . Recently, with the Chandra 
LETGS-HRC, Nicastro et al. (2002) discovered X-ray resonant absorption features from 
warm/hot local gas along the line of sight towards PKS 2155-304, giving the first X-ray 
evidence for possible WHIM gas. These absorbers are near zero redshift and therefore most 
likely associated with the Galaxy or local group (see § 2). 

We report the first detection of an O VIII Lyct absorption line along the sightline towards 
PKS 2155-304 at cz ~ 16,600 km s _1 plausibly associated with an intervening, modestly 
overdense region in the IGM. The observations were made with the Chandra LETG-ACIS 
(Advanced CCD Imaging Spectrometer). Physical diagnostics show the gas that gives rise 
to this resonance line has the typical properties of the WHIM, indicating the detection of a 
significant fraction of the "missing baryons" . 

2. Data Analysis 

PKS 2155-304 at z ~ 0.116 (Bowyer, Brodie, Clarke, & Henry 1984) is one of the 
brightest soft X-ray sources (e.g. Urry & Mushotzky 1982 and references therein) and has 
been used as a Chandra calibration target. It was observed with the Chandra LETG-ACIS 
(see, e.g., Chandra Proposers' Observatory Guide at http://asc.harvard.edu/) on May 31, 
2000 , December 6, 2000 and November 31, 2001. We restrict all spectra to the range 2 — 42 
A to avoid complexities at longer wavelengths from the instrumental C-K edge. The events 
were summed in 0.025A bins, half of the instrumental line response function, which has 
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nearly constant FWHM ~ 0.05 A across its bandpass (for data processing procedures, see 
McClintock et al. 2001). We corrected the wavelength scale by calibrating the spectrum with 
two known features (atomic O I from the Galaxy and solid-state O I from the instruments). 
The three data sets were separately divided by the best fit Galactic absorbed power-law to 
remove the underlying broad-band continuum behavior. The remaining residuals were then 
fit with a five-order polynomial. Such a technique removes spectral features larger than ~ 8A 
(e.g. 5-10% calibration uncertainties) but preserves narrow line features. The data from the 
three observations were then combined as a weighted sum that accounts for the differences 
in exposures and effective area. 

The observation during November 2001 has the strongest soft flux ( respectively for May 
2000, December 2000 and November 2001, F ai _ 2 .4 = 1.35, 0.87 and 2.92 xl0~ 10 ergs cm-V 1 ). 
All continua are well described by a single power law (r ~ 2.31 ± 0.25,2.60 ± 0.20, and 
2.72 ±0.47, respectively) absorbed by a Galactic hydrogen column of N H = 1-36 x 10 20 cm~ 2 
(Lockman & Savage 1995) (errors are quoted at 90% confidence). 

After a blind search for any statistically significant absorption features, two absorption 
features with S/N > 4 were detected in the 2-42 A region of the LETGS spectral bandpass 
(Figure 1). These features were subsequently fit in ISIS (Interactive Spectral Interpretation 
System : http:/ /space. mit.edu/ASC/ISIS/) with Gaussian models (Table 1). The absorption 
feature at ~ 21.6 A was reported by Nicastro et al. (2002) in the LETGS-HRC archival 
data, which they attribute to O VII Hea resonance absorption from warm/hot local gas. We 
also show in Table 1 our results for two other features claimed by Nicastro et al. (2002). 
The feature attributed to O VIII Lya is marginally detected (3a) while O VII He/3 is not 
significant (lcr). We defer the discussion of these three features to a later paper, noting 
only that they are at cz ~ and therefore probably associated with the Galaxy or local 
group. We searched for the 20A feature in the LETG-HRC data of this source, but find that 
this region is complicated by instrumental feature (s) that we also found in the LETG-HRC 
observations of 3C 273. 



3. An absorption feature at ~ 20 A — evidence of the WHIM 

3.1. Line Identification 

We concentrate on the absorption feature which appears at 20.02 A (619 eV). This falls 
at about the same energy where Canizares & Kruper (1984) detected a broad (~ 50-100 eV) 
absorption feature from an observation with the Einstein Observatory Objective Grating 
Spectrometer. A similar broad feature was detected by Madejski et al. (1991). However, the 
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line we detect is unresolved and has a width of < 0.039A (or 1.2 eV) at 90% confidence. The 
Poisson significance of that this feature is due to a chance fluctuation is 2.8 x 10 -6 , equivalent 
to 4.5(7 for Gaussian statistics. The probability of finding one such feature by chance in the 
1520 bins from 2-40 A is 0.35%. The 21.61 A feature has a similar confidence level, while 
the next most prominent feature has 20% probability of occuring by chance (these values 
were confirmed using monte carlo simulations). 

Considering cosmic abundances and oscillator strengths (e.g. Verner & Yakovlev 1995) 
for different ions, O VIII Lya is the only strong candidate line between 18 and 20 A , the 
measured wavelength de-redshifted to the source. One scenario which may reconcile O VIII 
Lya at 18.97A with the 20.02 A of the detected feature is that this is the signature of a 
high velocity ~ 18, 000 km s -1 wind/jet blueshifted from the BL Lac rest frame. In this case 
the narrowness of the line would be surprising, especially when compared to the very broad 
(pa 30,000 kms" 1 ) feature seen by Canizares & Kruper (1984) and Madejski et al. (1991), 
although we note that the relativistic jet in the galactic binary SS433 does have roughly 
constant velocity (Marshall, Canizares, & Schulz 2001). The implied mass outflow would be 
hundreds to thousands of M Q yr _1 , depending on assumptions about ionization conditions 
and beaming factors, as shown in Krolik, Kallman, Fabian, & Rees (1985). They also discuss 
some of the resulting difficulties with such an outflow. 

It is more plausible that the 20 A absorption is due to O VIII Lya in a known intervening 
system at cz pa 16, 734 km s _1 . With the Hubble Space Telescope (HST), Shull et al. (1998) 
discovered a cluster of low metallicity H I Lya clouds along the line-of-sight (LOS) towards 
PKS 2155-304, most of which have redshift between cz = 16, 100 km s _1 and 18, 500 km s _1 . 
Using 21 cm images from the Very Large Array ( VLA), they detected a small group of four 
H I galaxies offset by ~ 400 — 800 hfp 1 kpc 1 from the LOS, and suggested that the H I 
Lya clouds could arise from gas associated with the group. The velocity centroid of the four 
galaxies is (V ga i) = 16, 853 km s _1 , which within the uncertainties is identical to that of our 
absorption feature. Therefore the hot gas giving rise to the O VIII absorption can plausibly 
be associated with the same overdensity indicated by the H I Lya absorbers and the small 
galaxy group. 

Taking the absorption line to be O VIII Lya, the Doppler 6-parameter must be b = 
y/2a < 830 km s" 1 (90% confidence) since the line is unresolved. If b > 200 km s 1 , then the 
line is unsaturated and N(0 VIII) ~ 9.5 x 10 15 cm~ 2 . At b < 100 km s _1 the line is heavily 
saturated, but this can be ruled out by the lack of higher Lyman series lines in our spectrum. 



1 We use H = 70h 70 km s 1 Mpc 1 , and a standard ACDM model with Q m — 0.3 and SI a = 0.7 
throughout the paper. 
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Associating the absorber with the galaxy group detected by Shull et al. (1998), which has 
a radial velocity dispersion of (<T gal ) = 325 km s -1 , suggests that the unsaturated case is 
most appropriate, and we adopt that for the subsequent discussion. However, simulations 
assuming complete viralization of the hot gas predict b < 100 km s _1 for lines such as this 
one (Fang, Bryan & Canizares 2002), so saturation cannot be ruled out. 

3.2. Density and Temperature Constraints 

We can constrain the density of the absorbing gas, assuming it is associated with the 
intervening galaxy group. First, the upper limit on the line width of o ~ 580 km s _1 sets 
an upper limit to the path length of ~ Sh^Mpc, since otherwise the differential Hubble 
flow would excessively broaden the line. For N(0 VIII) ~ 9.5 x 10 15 cm~ 2 , this gives 
rib > (1-0 x 10~ 5 cm" 3 ) Zq I/q^I^ 1 where Z .i is the metallicity in units of 0.1 solar abundance, 
/o.5 is the ionization fraction in units of 0.5 and Z 8 is the path length in units of Sh^Mpc. We 
adopt a metallicity of 0.1 Z (Anders & Grevesse 1989), which is the mean value predicted 
from hydrodynamic simulations for overdensities of ~ 30, with a dispersion of a factor of 
two (Cen & Ostriker 1999b). A more reasonable estimate of the path length comes from 
the mean projected separation of ~1 Mpc for the galaxies in the group (Shull et al. 1998), 
which gives rib ~ 7.5 x 10~ 5 cm~ 3 Z^lf^l. In fact, this could be taken as a plausible upper 
limit, since the WHIM associated with this overdensity could well be more extended than 
the galaxies, as are the H I Lya clouds. This implies a range of baryon overdensity over 
the cosmic mean (n b ) = 2.14 x 10~ 7 cm~ 3 (Tytler, O'Meara, Suzuki, & Lubin 2000) of 5 b ~ 
50 - 350. Interestingly, Shull et al. (1998) estimate an overdensity for the galaxy group of 

Sgal ~100. 

In the case of pure collisional ionization, temperature is the only parameter of impor- 
tance over a wide range of density so long as the gas is optically thin. The O VIII ionization 
fraction peaks at 0.5, and exceeds 0.1 for temperatures T ~ 2 — 5 x 10 6 K. Using CLOUDY 
(Ferland et al. 1998) we find that photoionzation by the cosmic UV/X-ray background is 
not important for rib > 10~ 5 cm -3 . 

We can compare this estimate for T to what we might expect from the dynamics of the 
group of galaxies, if it forms a bound system. Then its virial mass is M vir « 3 x lO^hfo 1 M . 
Scaling the mass-temperature relationship down from richer clusters of galaxies (Bryan & 
Norman 1998), gives a temperature of ~ 9 x 10 6 K. Similarly, using the relation between 
velocity dispersion and temperature obtained from fitting a sample of galaxy groups gives a 
value of 1 - 2 x 10 7 K (Xue & Wu 2000). These values are somewhat higher than what we find 
from the ionization fraction. This could suggest that the system is not bound and virialized. 
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However, Mulchaey, Mushotzky, Burstein, & Davis (1996) suggested that the intragroup gas 
in spiral-rich groups, such as this one, may be cooler than that of groups and clusters of 
early-type galaxies. They suggested that these cooler systems may imprint absorption lines 
from highly-ionized species (e.g. O VI) on quasar spectra, which may be what we observe. 

CLOUDY calculations of the column density ratios between other ions and O VIII give 
another tight constraint on gas temperature. Based on the non-detection of other ion species, 
we calculate the 4<r upper limits of the column densities of H- and He-like C, N, Ne and 
He-like oxygen relative to O VIII. N(0 VII)/N(0 VIII) provide the strongest constraints on 
gas temperature. At roughly T > 10 6 4 K, all the column density ratios satisfy the observed 
constraints. 

Adopting the value of T = 4 — 5 x 10 6 K, we can compare the cooling timescale for the 
gas to the Hubble time tu- The cooling time scale t coo \ ~ (3/2)(kT/ribA(T)), here A(T) is 
the cooling function (Sutherland & Dopita 1993), which is a function of temperature and 
metallicity (in the temperature range of interest, line emission from iron is the most effective 
coolant). For 0.1Z Q , the baryon density should be n b < 2 x 1CT 4 (t H /t coo i) cm~ 3 so that 
the gas remains hot. This density limit, however, would be lower by a factor of ~ 3 if the 
abundance is solar. Alternatively, we can use the cooling argument to set an upper limit on 
the metallicity. Adopting a baryon density of 7.5 x 1CT 5 cm~ 3 and requiring t coo i > tn, the 
cooling function A must be less than 2 x 10~ 23 ergs cm 3 s _1 , which implies a metallicity of 
< O.3Z . 

Following Lanzetta et al. (1991) and Tripp & Savage (2000) we estimate the baryonic 
content f2b(0 VIII) that can be probed by our O VIII resonance absorption line, expressed 
in units of (rib). We assume conservative upper limits of Z < 0.5Z Q and O VIII ionization 
fraction of / < 0.5. Given the path length of Az < 0.116, we estimate £\(0 VIII) > 
O.OOShfo 1 . This is about 10% of the total baryon fraction, or about 30-40% of the WHIM 
gas, if the WHIM gas contains about 30-40% of total baryonic matter (Dave et al. 2001). 
This baryon fraction is consistent with the prediction from Perna & Loeb (1998) based on 
a simple analytic model; however, these numbers are very rough estimates, due to the large 
uncertainties in the abundance and ionization fraction. 



4. Discussion & Summary 

Our detection of O VIII Lya from an intervening overdense region in the IGM at cz ~ 
16,600 km s -1 , together with the detection by Nicastro et al. (2002) of absorption from 
systems at cz ~ associated with the Galaxy or local group, have begun to reveal the 
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much-anticipated warm-hot component of intergalactic matter (Hellsten, Gnedin, & Miralda- 
Escude 1998; Cen & Ostriker 1999a; Dave et al. 2001; Fang, Bryan & Canizares 2002). 

The O VIII Lya feature we observe is plausibly associated with the same cosmic over- 
density that gives rise to the H I Lya clouds and the four spiral galaxies found by Shull et 
al. (1998). Three of the galaxies have radial velocities within ~ 100 km s _1 of one another, 
while the fourth is within ~ 400 km s -1 . In contrast, the velocities of the H I Lya clouds 
span ~ 1800 km s _1 , with two falling very close to the velocities of the closest projected 
galaxies. The galaxies represent an overdensity of 5 ga i ~ 100 relative to the mean. Our 
estimates for the overdensity of the hot absorbing gas are 5b ~50-350. The H I Lya clouds 
have very low metallicity, indicating that they are primarily composed of primordial mate- 
rial, whereas the hot gas should have an oxygen abundance of ~ 0.1 solar to give sensible 
values for the column density. It is unlikely that the same gas is responsible for both H I 
and O VIII absorption. Following Shull et al. (1998), we identify this as a region of modest 
overdensity (the galaxy group) connected to filamentary or sheet-like structures demarcated 
by the H I Lya clouds, such as those seen in numerous cosmological simulations. 

We note that the detection of a system such as this is roughly consistent with previous 
simulations. Based on the models of Perna & Loeb (1998) and Fang, Bryan & Canizares 
(2002), we estimate the cumulative absorption line number along a random sightline (Fig. 2). 
We adopt the numerical relation between baryon overdensity and metallicity from Cen & 
Ostriker (1999b) and assume collisional ionization. This observation is ~ 1.5a above the 
predicted value at the measured column density. According to Fig. 3, the probability of 
detecting such an O VIII Lya line or stronger to z ~ 0.1 is ~ 10%. The consistency, if borne 
out by further studies, suggests that the simulations are valid descriptors of the warm/hot 
component of the IGM, at least for overdensities of ~100. Using those simulations, we 
estimate that we can probe about 10% of total baryons and about 30-40% of the WHIM gas 
in the local universe. 

Other observations of bright AGN with the Chandra and XMM-Newton grating spec- 
trometers are likely to reveal additional WHIM absorption features, and may turn the present 
few samplings into a real forest of high ionization Lya or Hea lines. However, the sensitivity 
of these instruments permit us only to probe the high density tail of the distribution (e.g. 
Fang & Canizares 2000). It will be left to future missions like Constellation-X and XEUS, 
assuming they have both high throughput and spectral resolving powers of 1000 or more at 
energies 0.1-1 keV, to fully reveal the X-ray forest. 

We are grateful to F. Nicastro and his colleagues for sharing an advance copy of their 
paper. We thank the other members of the MIT/CXC team for its support. We also thank 
the referee, J. M. Shull, for helpful comments. This work is supported in part by contracts 
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Table 1: Fitting parameters of the X-ray absorption Lines 



O vin Lya O vm Lya O vii Hea O vn Hc/3 



"rest 


18.9689 


18.9689 


21.6019 


18.6288 


^obs 


18 95+ - 02 


20 02+ 015 


21.611HJ 


18.63lg;g 


cz(km s _1 ) 


-300±l\l 


16624 ± 237 


112tl« 




Line Width 6 




< 0.039 


< 0.027 




Line Flux c 




4.8t?j 


tc rc+3.0 
°-°-1.7 


1 1+ 1 ' 8 
i - 1 -0.6 


EW (mA) 


o 0+3.5 
°- S -6.9 




15.6l|| 


o q+5.5 
°-* > -1.8 


SNR d 


3.0 


4.5 


4.6 


1.1 



a. Rest-frame wavelengths in units of A (Verner & 
Yakovlev 1995). 

b. 90% upper limit of the line width a, in units of A. 

c. Absorbed line flux in units of 10~ 5 photons cm _2 s _1 . 

d. In units of equivalent sigma of a Gaussian distribution 
with the same confidence level as the Poisson significance 
gives. 
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Fig. I.— The Chandra LETG-ACIS spectrum of PKS 2155-304 between 16 and 23 A. (a): 
The red solid line is the fitted continuum (Galactic absorbed power law plus the polynomial), 
the averaged la error bar plotted on the right is based on statistics only; (b) the signal-to- 
noise ratio; (c) enlarged portion of the spectrum between 18 and 22 A, the red solid line is 
the fitted continuum plus three Gaussian models, the ion species are labeled in blue. 

Fig. 2. — Cumulative O VIII absorption lines per unit redshift vs. column density. Three 
data points are: PKS 2155 (this paper), H 1821 (Fang et al. 2002), S5 0836 (Fang, Marshall, 
Bryan, & Canizares 2001). The errors are estimated at 2a level (Gehrels 1986). For non- 
detection cases, 2o" upper limits are given. 
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